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Following the ideas behind the Feynman approach, a variational wave 

function is proposed for the Frohlich model. It is shown that it provides, 

for any value of the electron-phonon coupling constant, an estimate of the 

polaron ground state energy better than the Feynman method based on path 

integrals. The mean number of phonons, the average electronic kinetic and 

interaction energies, the ground state spectral weight and the electron-lattice 

correlation function are calculated and successfully compared with the best 

available results. 
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I. INTRODUCTION 



In recent years a large amount of experimental results has pointed out that the electron- 
phonon (e-ph) interaction plays a significant role in determining the electronic and mag- 
netic properties of new materials as the high Tc superconductors and the colossal magneto- 
resistance manganites.^ The experimental data have given rise to a renewed interest in 
models of the e-ph coupled system. In this paper we investigate the polaronic features of 
the Frohlich model within a variational approach.^ Here the picture is the following. When 
an electron in the conduction band of a polar crystal moves through the crystal, its Coulomb 
field produces in its neighborhood an ionic polarization that will infiuence the electron mo- 
tion. Then the particle must carry this polarization with it during its motion through the 
crystal. The quasi-particle formed by the electron and the induced polarization charge is 
called polaron. Within the Prohlich model: 1) the optical modes have the same frequency; 
2) the dielectric is treated as a continuum medium; 3) in the undistorted lattice the electron 
moves as a free particle with a quadratic dispersion relation (effective band mass approxi- 
mation) . 

The problem of finding the ground state energy of the Frohlich Hamiltonian attracted 
the interest of a lot of researchers mainly in the period 1950-1955. Numerous mathemati- 
cal techniques have been used to solve this problem: from the perturbation theory in the 
weak coupling regime^ to the strong coupling theory,"^ from the linked cluster theory^ to 
variational^ and Monte Carlo approaches. ^'^ The weak coupling regime is well described 
within the Lee, Low and Pines (LLP) approach.^ Here, after the dependence of the Hamil- 
tonian on the electron coordinates has been eliminated, an upper bound for the polaron 
ground state energy is obtained by using a variational wave function which is based on the 
physical assumption that successive virtual phonons are emitted independently. In the op- 
posite regime, when the e-ph interaction is very strong, a good description of the polaron 
features has been obtained by Landau and Pekar.^'^ Their theory, based on a variational 
calculation, stems from the idea that, for very large values of the e-ph coupling constant. 
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the electron can follow adiabatically the quantum zero-point fluctuations of the polarization 
field. In their first papers the electron is localized with a Gaussian wave function. Next, 
the method has been improved by Hohler^^ by constructing an eigenstate of the total wave 
number by superposing Landau-Pekar states localized at different points of the lattice. In 
any case the validity of LLP and Hohler approaches is restricted, respectively, to weak and 
strong e-ph coupling regimes. 

An excellent approximation, that is accurate at all couplings, has been introduced by 
Feynman.^^ His approach provides a variational estimate of the electron self-energy based 
on the path integrals and the Feynman- Jensen inequality. After the phonon variables have 
been eliminated exactly, Feynman introduces a model Hamiltonian which describes approx- 
imatively the interaction of the electron with the lattice. This Hamiltonian is that of an 
electron coupled to another particle with a harmonic oscillator coupling. The trial action for 
the system is obtained by eliminating the coordinates of the fictitious particle simulating the 
phonon degrees of freedom. The mass M of the fictitious particle and the spring constant 
are the two variational parameters within the Feynman approach. The Monte Carlo study^'^ 
of the Frohlich model has demonstrated the remarkable accuracy of the Feynman method 
to the electron self-energy. 

In this paper we use a variational technique, within an Hamiltonian approach, to investi- 
gate the polaronic features of the Frohlich model. It is based on linear superposition of two 
translationally invariant wave functions that provide a very good description of the weak and 
strong e-ph coupling regimes. These wave functions are built assuming as starting points 
the LLP^ and Hohler^^ variational approaches. First, we improve these methods obtaining 
a better upper bound for the polaron ground state energy in the two asymptotic regimes of 
weak and strong e-ph interaction. Then, we use a linear superposition of these two wave 
functions. The comparison of our results with the Feynman^^ and Monte Carlo data^ shows 
that the proposed method provides an excellent description of the polaron ground state 
energy for any value of the e-ph coupling. Within our variational approach, the estimate of 
the electron self-energy turns out systematically lower than one of the Feynman method. In 
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particular, unlike the Feynman approach, the ground state energy shows the exact depen- 
dence on the e-ph coupling constant in the strong coupling regime. Next, we calculate the 
mean number of phonons present in the virtual phonon cloud surrounding the electron, the 
average electronic kinetic and interaction energies, the ground state spectral weight and the 
induced ionic polarization charge density. These quantities are successfully compared with 
Monte Carlo® and Feynman results. 

The proposed method has the advantage to exhibit, first to the author's knowledge, 
a wave function that gives the correct behavior in both weak and strong coupling limits 
and provides an interpolation between them with results at least accurate as those of the 
Feynman approach. 



II. THE MODEL 

The Frohlich modeP is described by the Hamiltonian: 

H^Hei + Hph + H,_ph = §^ + T.^oala^+ X^(M,e^«"'-%-+ h.c). (1) 

q q 

In Eq.(l) m is the band mass of the electron, hujo is the longitudinal optical phonon 
energy, r and p are the position and momentum operators of the electron, at represents 
the creation operator for phonons with wave number q and Mg indicates the e-ph matrix 
element. In the Frohlich model,^ Mg assumes the form: 



Mg^ihuJo^J-^ , (2) 



where a, dimensionless quantity, is the e-ph coupling constant, Rp = ^^2:;^ and V is the 
volume of the system. 
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III. THE STRONG COUPLING REGIME 



A. The adiabatic approximation 

When the value of a. is very large (a ^ 1) the electron can follow adiabatically the lattice 
polarization changes and it becomes self-trapped in the induced polarization field. The idea 
of Landau and Pekar,^° in the first works on polarons, is to use, as trial wave function for 
the e-ph coupled system, a product of normalized variational wave functions |(/?) and |/) 
depending, respectively, on the electron and phonon coordinates: 

l^) = k)l/)- (3) 
The expectation value of the Hamiltonian (1) on the state (3) gives: 

2 

= + (/I E [^^oat«<r+ + P;4] I/) (4) 



with 



The variational problem with respect to |/ > leads to the following lowest energy phonon 
state: 

|/>=e^n''"o"^ '^•"J|0> . (6) 

The minimization of the corresponding energy with respect to |(^) leads to a non-linear 
integro differential equation that has been solved numerically by Miyake.^^ The result for 
the polaron ground state energy in the strong coupling limit is: 

E = -0.108513a^/iu;o- (7) 
The Landau-Pekar^° Gaussian ansatz for \ip): 



after the minimization of the expectation value of the Hamiltonian (1) on this state, with 
respect to the variational parameter cu, provides an estimate of the ground state energy: 
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E = -—huo ^ -0.106103q;^;iu;o (9) 

OTT 

that is very close to the exact result (7). The best value for u; turns out: 

u = -^ujQ. (10) 

An excellent approximation for the true energy (7) is obtained by using a trial wave 
function similar to that one introduced by Pekar:^'^ 

\ipp) = Ne-'^'' [l + b (27r) + c (11) 

with N normalization constant and b, c and 7 variational parameters. The minimization of 

{ipp\H\ipp) leads to: 

E ^ -0.10S507a'^nuJo- (12) 
This upper bound for the energy differs from the exact value less than 0, 01%. 

B. Path integral method versus Hamiltonian approach 

At this point we recall the result of the Feynman^^ variational calculation when the 
approximating action is represented by a fixed harmonic binding potential: 

hcuo, a 00. (13) 



E 



-— -31og2 
37r ^ 



It is given by the sum of two terms. The first one corresponds to use a Gaussian wave 
function in the Landau and Pekar's method (Eq.(8)). The last one does not depend on 
the e-ph coupling constant a. The origin of this contribution in Feynman's expansion of 
the polaron energy has been discussed by AUcock^^ by using the perturbation theory in the 
strong coupling limit. Our first aim is to put this result on variational basis. This will allow 
us to characterize the terms that one has to introduce in the trial wave function to improve 
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the Landau and Pekar's ansatz. To this aim, starting from Eq.(3) and Eq. (6), we apply 
the following unitary transformation: 



(14) 



with 



- — ag — h.c. 



(15) 



The transformed Hamiltonian assumes the form: 



Hi — Ho + Hi 



(16) 



with 



T) / \ CX-' 

^0 = 1^ + E ^^oaJ-«,^ - E {M.e"^''' - a^) ^ + h.c. 



huJn 



-E 



^ fKjJo 



(17) 



and 



(18) 



One recognizes immediately that the Landau-Pekar approach corresponds to use as trial 
wave function for Hi. 



(19) 



with \ipip) given by Eq.(8) and ag = Mg{ipip\e^^'^\(pip) . In other words, in this approach, one 
approximates the lowest energy state of Hq with a Gaussian wave function containing the 
variational parameter cu, that represents the characteristic oscillation of the electron in the 
induced lattice polarization. The next order term is obtained assuming Hj as perturbation 
and approximating the eigenstates of Hq with those of an harmonic oscillator. At the first 
order of the perturbation theory the wave function is: 



(20) 



where 

Mr, t) = ^e^^^-r\4^i^'-^) - ^ . (21) 

This expression has been got using the generating function of the Hermite polynomials. 
Finally we note that can be obtained from 

|V') = e-^^|v^,p)|0) (22) 

with 

82^- 1 y. [hq{r, t)aq - h.c] dt , (23) 

by expanding the exponential e~^^ and truncating the expansion at the first order. Taking 
into account also the unitary transformation in Eq.(14), the previous considerations lead us 
to assume as trial wave function for the Prohlich Hamiltonian in the strong coupling limit: 

\^F)^e-'-\cpi,)\0), (24) 



with 



Vf«eV2^(*^-i)a„— /i.e. 



dt . (25) 



We have indicated this coherent state with " F" since it is easy to show that the expectation 
value of the Prohhch Hamiltonian on the state (25) gives: 

i.e. the Feynman result when the approximating action is represented by a fixed harmonic 
binding potential. In Eq.(26) T[x) is the Gamma function. In particular the minimization 
of E with respect to the variational parameter uj and the asymptotic expansion for a —> 00 
restore the Eq.(13). Then, the order beyond the Landau and Pekar's theory is due to the 
lattice fluctuations and to the consequent change in the electron wave function. 
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C. Improvements of the Feynman result 



The next step is to try to improve the Feynman result. To this aim, we note that is 
possible to obtain an excellent approximation of the polaron ground state energy in Eq.(26) 
substituting in Eq.(25) Sp with: 



where 



huj Jo 



(27) 



(28) 



and 



Un 



flUJ 



f 

J a 



J g 2 2muj 



(29) 



Here a and r] are two variational parameters. In other words, we obtain the main contribu- 
tion to the Feynman estimate of the electron self-energy approximating Sp as sum of two 
terms: the first one stems from the observation that the electron moves very fast in the 
induced potential well; the second one takes into account the lattice fluctuations and the 
possibility that they can follow instantaneously the electron motion. In order to improve 
the Feynman result, the Pekar's approach (Eq.(ll)) and the previous analysis suggest us to 
try the following ansatz: 



(30) 



with |(/9p) given by Eq.(ll), rj variational parameter and /^-and s^- functions to be determined 
by minimizing the expectation value of the Frohlich Hamiltonian on this state. This last 
quantity turns out: 



2^2 



2m 



+ E 



huJo + J {r^s^rt + h.c) - (M^st + M.r^t + h.c) (31) 



with 

r,-=((^,|e^'?-(i-'')|(^^). (32) 

Making {tplHlip) stationary with respect to arbitrary variations of the functions and 
Sg-, we obtain two, easily solvable, algebraic equations. The minimization and the asymptotic 
expansion of the ground state energy, for a — > oo, provide: 



E 



-0.108507a2 - 1. 



fiujQ. (33) 



The electron self-energy shows the exact dependence on in the strong coupling regime 
together with a good estimate of the e-ph coupling constant independent contribution due to 
the lattice fluctuations. This allows to obtain, for a > 8.7, an upper bound for the polaron 
ground state energy better than the Feynman approach when the approximating action is 
represented by a fixed harmonic binding potential (Eq.(26)). On the other hand, both these 
methods give the same result for a <6, i.e. E — —ahcuo- However, both the methods show 
a discontinuity in the transition from the weak to strong coupling regime. 

To overcome this difficulty one has to take into account the translational invariance. We 
construct an eigenstate of the total wave number by taking a superposition of the localized 
states (30) centered on any point of the lattice in the same manner in which one constructs 
a Bloch wave function from a linear combination of atomic orbitals: 

l^(.c)) = JiP{f-R)d^R. (34) 

The minimization, with respect to the variational parameters, of the expectation value of 
the Prohlich Hamiltonian on this state, that accounts for the translationally symmetry, and 
the asymptotic expansion for a — > cxo provide: 



E 



-0.108507q;^ - 2.67 



hcuQ. (35) 



This upper bound is lower than the variational Feynman estimate which for large values of 
a assumes the form:^^ 



E 



a' , 3 
---31og2-- 



hcuQ. (36) 



10 



IV. THE WEAK COUPLING REGIME 



When the value of a is very small the lattice follows adiabatically the electron. A good 
physical description of the polaron features in this regime is provided by the LLP approach.^ 
The starting point is the observation that the total momentum operator: 

Pt=P + Y^ hqaia^ (37) 
q 

is a motion constant, i.e. it commutes with the Hamiltonian. The conservation law of the 
total momentum is taken into account through the unitary transformation: 

^^^<Q-E.^«W)-^ (38) 

where HQ is the eigenvalue of Pt- In this paper we are interested in the ground state 

— * 

properties of the e-ph coupled system, so that we will restrict ourselves to the case Q — 0. 
The transformed Hamiltonian does not contain the electron variables and it is given by: 

(39) 

The LLP wave function is: 

\^)^^lf'i'^^-^-<^)\^)^ (40) 

where |0) is the phonon vacuum state and — M^j (huo + In other words 1-0) is 

the lowest energy state of the first two terms of the transformed Hamiltonian Hi. The use 
of this wave function is based on the physical assumption that, when the e-ph interaction 
is weak, there is not correlation among the emission of successive virtual phonons by the 
electron. This assumption restricts the validity of this approach to the regime characterized 
by small values of a. The ground state energy turns out E = —ahuQ. In other words, this 
method puts the results of the perturbation theory on variational basis. 

To improve the LLP approximation,^ one has to introduce in the trial wave function a 
better description of the recoil effect of the electron, effect present only on average in LLP 
approach. This can be done using the following ansatz: 
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c)) = e 



|o) + E ^9i,9-i4i4Jo) 



li at 



(41) 



that takes into account the correlation between the virtual emission of pairs of phonons.^^ 
In this paper we will choose: 

and 

Here 7, 6 and e are three variational parameters that have to be determined by minimizing 
the expectation value of the Hamiltonian (1) on the state (41). This procedure provides as 
upper bound for the polaron ground state energy at small values of a: 

E = -ahuQ - QM2?,a^rwjQ, a ^ , (44) 

i.e. the same result, at this order, of the Feynman approach. We stress that, at the 
order, the result for the electron self-energy is: 

E = -ahcuo - OM59a^nuJo (45) 

as found by Hohlcr and MuUensicfcn,^'' Larsen^® and Roscler.^^ 



V. ALL COUPLINGS 

A careful inspection of the wave function (34) shows that is able to interpolate between 
strong and weak coupling regimes. On the other hand, for small values of a a better 
description of the polaron ground state features is provided by the wave function (41). 
Moreover, in the weak and intermediate e-ph coupling, a < 7, these two solutions are not 
orthogonal and have non zero off diagonal matrix elements. This suggests that the lowest 
state of the system is made of a mixture of the two wave functions that give an accurate 
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description of weak and strong e-ph coupling regimes. Then the idea is to use a variational 
method to determine the ground state energy of the Hamiltonian (1) by considering as trial 
state a linear superposition of the two previously discussed wave functions: 



where 

' \wc)) > _ l'^(sc)) 



(46) 



i^(~))= ,,;7 . . i^(..)>= . (47) 



and S is the overlap factor: 



In Eq.(46) A and B are two additional variational parameters that provide the relative 
weight of the two solutions in the ground state of the system. In this paper we perform 
the minimization procedure in two steps. First, the expectation values of the Prohlich 
Hamiltonian on the two trial wave functions in Eq.(41) and Eq.(34) are minimized and the 
variational parameters are determined. Then, the minimization procedure discussed in the 
present section is carried out. This way to proceed simplifies significantly the computational 
effort and makes all described calculations accessible on a personal computer. An approach, 
similar to that one described in this section, has been successfully used for the Holstein 
model. 

The procedure of minimization of the quantity {ip\H\ip) with respect to A and B gives 
for the polaron ground state energy 



E ^ V _ ^ L (43) 



and for the ratio of the two parameters A and B 



Ec-ES 



B E — E(^wc) 

Here = {■ip(wc)\H\'^{wc)) , ^(^c) = (^(sc)l^l^(sc))> = (^(«;c) + ^(.c)) /2 and E, 

iwc)m(sc)) + h.c.) /2. 
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(50) 



VI. NUMERICAL RESULTS 



In Fig.l we plot the polaron ground state energy, obtained within our approach, as 
a function of the e-ph couphng constant a. The data are compared with the results of 
the variational treatments due to Lee, Low and Pines, ^ Pekar,^° Feynman^^ and with the 
energies calculated within a diagrammatic Quantum Monte-Carlo method/ As it is clear 
from the plots, our variational proposal recovers the asymptotic result of the Feynman 
approach in the weak coupling regime, improves the Feynman's data particularly in the 
opposite regime, characterized by values of the e-ph coupling constant a ^ 1, and it is 
in very good agreement with the best available results in literature, obtained with the 
Quantum Monte Carlo calculation/ This agreement indicates that the true ground state 
wave function is very close to a superposition of the above introduced functions, that provide 
a very good description of the two asymptotic regimes. Within our approach we have also 
calculated the mean number of phonons present in the virtual phonon cloud surrounding 
the electron, A^, the average electronic kinetic and interaction energies, K and /. These 
quantities are reported, respectively, in Fig. 2, Fig. 3 and Fig.4 where are compared with the 
same properties obtained in the Feynman's variational treatment based on path integrals. 
As for the ground state energy, our variational ansatz is able to recover all the expected 
behaviors. For small values of a, A^ ^ q;/2, K aTiujQ/2, I —2ahujQ as predicted by 
the LLP approach^ and the weak coupling perturbation theory. In the opposite regime 
the electron follows adiabatically the lattice polarization changes. The values N — 

2 A '2 

K = ^huo, I = —^huiQ obtained within the Landau and Pekar's variational treatment 
(see Eq.(8)), based on the electron self-trapping with a Gaussian wave function, represent 
very accurate estimates of these quantities when they are calculated within the Feynman's 
approach. On the other hand, the values N = 2 - 0.108507q;^ K = O.lOSSOTa^^a;^ / = 
— 4-0.108507Q;^/ia;o obtained within the Pekar's variational treatment (see Eq.(ll)) represent 
very good approximations for the same quantities calculated within our approach. Then the 
variational Feynman's and our methods differ mainly in the strong coupling regime as it 
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turns out from the plots in Fig.l, Fig.2, Fig.3 and Fig.4. We stress that, in this range of 
values of a, our approach provides a better estimate of the polaron ground state energy than 
the Feynman's method. 

Another physical quantity of interest is p(r), i.e. the average ionic polarization charge 
density induced at a distance r by the electron. This quantity is related to the static 
correlation function between the electron position fl = and the oscillator displacement at 



r: 



= - (^) im = 0)1 Y.iM,e'^"'Q\+ h.c.mre = 0)). (51) 
It easy to show, analytically, that the exact sum rule for the total induced charge:^° 

/ p{r)d'r = e (— - i) (52) 

J V Coo ^0 / 

is satisfied within our variational approach. Figure 5 shows p(r)/ / p{r)d^r as a function of 
r for different values of the e-ph matrix element a corresponding to weak, intermediate and 
strong coupling regimes. Our data are compared with results obtained within the Feynman's 
method^^ and a path integral Monte Carlo scheme.^° If the e-ph couphng is weak, the lattice 
deformation is not able to trap the charge carrier. The extension of the polaron is large 



compared with the characteristic length y 2mujo ■ situation is different in the opposite 
regime where the lattice deformation is localized around the electron. In any case also this 
correlation function, evaluated within our approach, is in agreement with the best data 
available in literature. 

Finally Figure 6 shows the ground state spectral weight: 

where |0) is the electronic vacuum state containing no phonons and ct is the electron creator 
operator in the momentum space. Z represents the renormalization coefficient of the one- 
electron Green's function and gives the fraction of the bare electron state in the polaron 
trial wave function. This quantity is compared with that one obtained in the diagrammatic 
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Quantum Monte Carlo method/ The result of the weak couphng perturbation theory is also 
indicated: Z = l — a/2. For small values of a the main part of the spectral weight is located 
at energies that correspond approximatively to the bare electronic levels. Increasing the 
e-ph interaction, the spectral weight decreases very fast and becomes practically zero in the 
strong couphng regime. Here the most part of the spectral weight is located at excited states. 
The diagrammatic Quantum Monte Carlo study^ of the Frohlich polaron has pointed out 
that there is no stable excited states in the energy gap between the ground state energy and 
the continuum. There are, instead, several many phonon unstable states at fixed energies: 
Ef — Eq ~ 1, 3.5 and S.bhujQ. These results seem to be contrary to the data about the 
optical absorption of large polarons,^^ which show, for large values of a, the presence of a 
very narrow peak corresponding to the electronic transitions from the ground state to the 
first relaxed excited state (RES). The nature of the excited states and the optical absorption 
of polarons in the Frohlich model require further study which is beyond the scope of this 
paper. 

In conclusion, in this paper, a variational approach has been developed to investigate 
the features of the Frohlich model. It has been shown that a hnear superposition of two 
wave functions, that describe the two asymptotic regimes of weak and strong e-ph coupling, 
provides an estimate of the polaron ground state energy which is in very good agreement 
with the best available results for any value of the e-ph matrix element. All the evaluated 
ground state properties show that the crossover between the two asymptotic regimes is very 
smooth. On the other hand the transfer of spectral weight from the polaron ground state to 
the higher energy bands turns out very fast. We stress that, to the best of our knowledge, it 
is the first time that a variational wave function, able to interpolate between the weak and 
strong e-ph coupling regimes, at least carefully as the Feynman method, is exhibited for 
the Frohlich model. 
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FIGURE CAPTIONS 



Fig.l. (a) The polaron ground state energy, E, is reported as function of a in units of 
hcuQ. The data (sohd hne), obtained within the approach discussed in this paper, are 
compared with the results (diamonds) of the Feynman approach, Ep, and the results 
(stars) of the diagrammatic Quantum Monte-Carlo method, Emc, kindly provided by 
A.S. Mishchenko. (b) The differences: E — Ep (diamonds) and E — Emc (stars) are 
reported as function of a. 

Fig.2. (a) The mean number of phonons, N, is plotted as function of a. The data, 
obtained within the approach discussed in this paper (solid line) , are compared with the 
results of the Feynman approach, Np (diamond), and the results of the diagrammatic 
Quantum Monte-Carlo method, Nmc (stars), extracted from Fig. 8 of ref.7. (b) The 
differences: Np — N (diamonds) and Nmc — N (stars) are reported as function of a. 
The error bars are due to uncertainty in the procedure used to extract the numerical 
values from Fig. 8. 

Fig. 3. (a) The average electronic kinetic energy, K, is plotted as function of a in units 
of Hcuq. The data (solid line), obtained within the approach discussed in this paper, 
are compared with the results (diamonds) of the Feynman approach, Kp. (b) The 
difference: Kp — K (stars) is reported as function of a. 

Fig. 4. (a) The average electronic interaction energy, /, is plotted as function of a in units 
of fiujQ. The data (solid line), obtained within the approach discussed in this paper, 
are compared with the results (diamonds) of the Feynman approach, Ip. (b) The 
difference: Ip — I (stars) is reported as function of a. 

Fig. 5. The average normalized ionic polarization charge density, induced at a distance r 
by the electron, is reported for three different values of a: (a) a = 1; (b) a = 6; (c) 
a — 12. The data (solid line), obtained within the approach discussed in this paper. 
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are compared with the results (dashed line) of the Feynman approach and the results 
(dotted line) of the Monte-Carlo method, kindly provided by S. Ciuchi. The distance 
r is measured in units of Rp. 

Fig.6. The ground state spectral weight, Z, is plotted as function of a. The data (solid 
line), obtained within the approach discussed in this paper, are compared with the 
results (stars) of of the diagrammatic Quantum Monte-Carlo method. The result of 
the weak coupling perturbation theory (dashed line) is also indicated. 
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